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Abstract 

Generalizing a result for the binary lens, similar alternative expressions 
are also given for the Jacobian determinant for a gravitational lens con- 
sisting of an arbitrary number of discrete lensing centres, with arbitrary 
masses and locations in a plane. The results for a binary lens are recov- 
ered, agreeing with some already known results, and correcting errors in 
others. 

1 Introduction 

Since the earliest studies of Einstein (1915,1916) on the deflection of light by 
gravitational fields in general relativity, a vast research and associated publica- 
tion record in gravitational lensing, both observational and theoretical, has blos- 
somed. The simple single mass point lens model is already the main paradigm 
for the understanding of the many 10^) lensing events detected by the var- 
ious monitoring and follow-up programmes under way. Binary lenses are also 
the focus of much effort, not only for their own intrinsic interest, but also for 
their potential as a tool, particularly in the detection of planetary systems. 
More sophisticated lens models, including many-body systems, are of interest 
too because of their possible galactic and cosmological applications. 



2 The lens equation(s) 

In an orthogonal rectangular co-ordinate system (x,y,z), with the origin at the 
observer, and the z-axis as optical axis, let a source be located in a (source) 
plane S at (angular size) distance Dg, and let there be n masses located in a 
(lens) plane L at distance Z^l- Rays from a source point at angular co-ordinates 
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{xs, Us), on crossing the lens plane are so deflected as to produce an image point 
(of which there may be more than one) at angular co-ordinates {xi,yi). The 
lensing conditions for small angles, see for example Paczynski (1996), are given 
by: 

Xi-Xs = a^{Ds - £'l)/-Ds; yi-ys = (Xy{Ds - -Dl)/-Ds 

where and ay are the angular deflections in the directions of the x and y 
axes respectively. In a two-dimensional vector notation, with r — (x, y) and 
<5 — (oix, Oiy) the above conditions may be written as the single equation ; 

For a spherically symmetric mass M the Einstein (1916) gravitational deflection 
is a = AGM/c^Rq, where Rq is the radial distance of closest approach or impact 
parameter. In the (rectilinear) ray optics approximation, the image position is 
deemed to be where the ray crosses the lens plane. Thus for each lensing mass 
Ml located at n (Z = l,n), setting fu = In-fil, ni = \fii\, and Rq = DlTu, the 
contribution to the deflection is di = ao{ri — ri)/r'ji where ag = ^GMi/ (?Di^. 
The lens equation for the n-body system then takes the form: 

n n 

ri-fs = [{Ds - £'l)/-Ds] 51 = 51 ~ 

1=1 1=1 

with mi = [{Ds - Di,)/DsDi}AGMi/c^ = {Mi/M)rl, where M = ^2^1 while 
rE = {4:GMc~^{Ds — Di,)/DsDi]^/^ is the Einstein iim(? angular radius for mass 
M. Without any loss of generality one may set m = ^ m; = 1, thus expressing 
all the angular measures in units of te. The lens equation may alternatively be 
written (Witt, 1990) in complex variable notation with z = x + iy and z = x — iy 
as : 

n n 

Zi - Zs = ^mi{zi - zi)/\zi - zi\^ = ^mi/{zi - zi) 
1=1 1=1 

3 Image amplification 

Since the intensity of radiation is invariant along a ray the brightness of the 
lensed source relative to that of the unlensed source (the magnification or am- 
plification) is the ratio of the areas or solid angles subtended at the observer, 
given hy A = J~^, where J is the Jacobian (determinant) of the transformation 
mapping fi into fg : 

J = \dfs/dfi\ = \d{xs,ys)/d{xi,yi)\ = \Sj;y - a^y\ 

with (ixx — Oiyyj o,xy — dyx and : 

axx = 53 mipY'^{r]f - a^^ = 2 mpi^^m 
I I 
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where = Xi ~ xi, rji = m - yi, and pi = m . Then : 

J = (1 - a^^){l + a^^) -aly = l- al^-aly = l-K 

I I 
Writing K = Ki + K2, expanding each term, and separating the 'diagonal' and 
'non-diagonal' parts of each, gives Ki = Kn + K12 and K2 = K21 + K22 where 

K^^ = EKpr'('?f - )]'; ^21 = Y.'y^mipi'^imf 

I I 

= 2 E E m^rnup-'pl\ri'i - eM - 

3^k k 

K22 = S^^m.jnikPj'^p^'^S.jVj^kVk 

j^k k 

Adding the 'diagonal' terms Ku and K21 gives : 

K3 = E "'-'^r^]^ ^ 2 E E 'f^j'^kPj^Pk^ = + 

I j^k k 

where K'32 may also be written : 

j^k k 

It then follows that : 

K,2 + i^32 = -4 E E ^o^kP-'p-AeA + ekv]] 

j^k k 

whereupon, adding the remaining 'non-diagonal' term if 22 gives K4 = K12 + 
K22 + K32 with : 

K4 = -4:^^mjmkPj'^p',:'^[£,jilk - £,kiljf 

i^k k 

SO that finally K = K^i + Kiov : 

K = ^mipl'^Y - '^Y^^^j^kP^'^ Pl^ijTlk - ^kVj? 

I i^k k 

The same result may naturally be obtained, perhaps more easily, from the 
complex variable formulation of the lens equation. There it may be shown 
(Witt, 1990) that : 

J = {dz,/dzif - \dzsidzi\^ = 1 - 1 E^^'Cr'r 
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where Q = Zi — zi and pi = \C.i\. Thus, with K = J — 1 : 

j k j k 

Expanding and separating the 'diagonal' and 'non-diagonal' parts gives K = 
K5 + Kq with : 

I I jjik k 

jjik k 

Putting Kr, = Kr,i + Kr,2, where ifsi = the second term may be 
rewritten : 

K52 = -2^^mjmkPj^p'^^[(j(jCk(k] = K32 

jjik k 

while Kq may be rewritten as : 

^6 = E E mjm,pj%^[(]ek + C]Cl] 

i^k k 

Now adding the 'non-diagonal' terms Kr,2 and Kq gives : 

-^''7 = E E ^j^kPj^Pk^iCjCk - Cj(k? = Ki 

jjik k 

since [CjCfe - CjCfe] = '^iiVj^k - ^jVk] thus giving the result K = K51 + K-j = 
K^i + K4, as before. It may be remarked here that [^j??fc — ^k'Hj] = ^ ^ife] so 
that : 

where the two-dimensional gradient Vj = d/dvi. 

4 The general binary lens 

In the case n = 2 the double summation reduces to a single term with j = l,k = 
2 giving : 

J = 1 - [mipf^ + m2P2 + 4mim2prV2 ^[Ci2]^ 
where C12 = [^i?72 — C2»?i] now reduces to : 

C12 = [{xiy2 - 2:22/1) + Xi{yi - y2) - yi{xi - X2)] 
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which may be simplified by an appropriate choice of co-ordinates. Taking the 
origin anywhere on the fine joining the two masses, then y\/x\ = y2/x2 and 
xiy2 - X2yi = gives C12 = [(2:2 - xi)yi — {y2 — yi)xi\ which is the same 
function of Xi and yi for given values of X2 — xi = 2xo and 2/2 — 2/1 = 2j/o, 
namely C12 = 2[a;o2/i — yoXi]- The further simplification of setting yo = 0, so 
that the x-axis lies along the line joining the two masses, with 7/1 = 7/2 = 0, 
gives C12 = 2xoyi. There still remains the choice of the value of Xi or X2, 
for which there are three 'natural' possibilities : Case (a) - Origin mid-way 
between the masses with xi = —xo,X2 = +xo; Case (b) - Origin at mass Mi 
with xi = 0,X2 = +2x0; Case (c) - Origin at mass M2 with xi = — 2a;o, X2 = 0. 

4.1 The critical curves 

The critical curves, or loci of infinite amplification, are given by the condition 
J = 0, which on rationalising takes the form D12 = with : 

D12 = p\p\J = pfpi - [mipl + m2pif + l&mim2xly'l 

in agreement with Schneider & Weiss (1986). Adopting plane polar co-ordinates 
= {ri,9i), then £'12 = dm cos'"(0i) where the coefficients Um depend on 
the choice of origin : 

Case (a) - fi = {ro,do) ■ 

«o = {rl + a;o)^ - {r^ + xl)'^{mi + 7712)^ + 16mim2rlxl 

ai = 4roXo{rl + xl){ml - ml) 
02 = -^rlxl[2{rl + xlf + (mi + 7722)^] 
as = 0; 04 = IGr^XQ 

In the symmetric case where mi = m2 = 1/2, ai = 0, the term in cos^j drops 

out, and the equation reduces to a quadratic equation in cos^ 9i, as given by 
Equation 9a of Schneider & Weiss (1986), with the solution : 

cos2 0,^[1 + 2{rl + xir + {1 + 8{4 + xt)y/']/{8rlxl) 

Case (b) - u = {ri,9i), with n = pi : 

ao = (r? + 4xlf{ri - m?) - 2r?(r? + 4a;^)mim2 

1 1 c 2 2 4 2 

-|-16ria;omim2 — rim2 
ai = 8riXo[rlmim2 - {rj + ^xl){rf - ml)] 
02 = 16rlxQ[rf — mf — mim2] 
Case (c) - ri = (r2, ^2), with r2 = P2 ■ 

ao = {rl + 4a;^)^(r| - mj) - 2rl{rl + Axl)mim2 
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, 1 (J 2 2 4 2 

+ lor2XQniini2 — 

ai = 8r2Xo[{rl + 4a;Q)(r2 — mg) — r2mim2] 

02 = 16r'^x'^[r2 — — mim2] 

from which it wih be evident that Case (b) and Case (c) are equivalent with 
respect to the interchange of the labels 1 and 2 together with the change of sign of 

Xq . Both therefore correct several errors in Equation 9b of Schneider and Weiss 
(1986), after allowing for their relative disposition of the binary masses. In the 
limit xo = the equations for all three cases degenerate to rf — (mi +m2)^ = 0, 

giving r-i = {mi + 777,2)^^^, so that the critical curve is the Einstein Ring for the 

combined masses of the binary. On the other hand in the limit xo = go, where 

either pi = 00 or p2 = 00, the condition J = reduces to pf — mf = 0,{l = 

1/2 

1,2), with the solution pi = rrii representing the now detached Einstein Rings 
around each mass. 
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